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Bending rigidity of a charged membrane or polyelectrolyte screened by monovalent counterions is 
known to be enhanced by electrostatic effects. We show that in the case of screening by multivalent 
counterions the electrostatic effects reduce the bending rigidity. This inversion of the sign of the 
electrostatic contribution is related to the formation of two-dimensional strongly correlated liquids 
(SCL) of counterions at the charged surface due to strong lateral repulsion between them. When 
a membrane or a polyelectrolyte is bent, SCL is compressed on one side and stretched on the 
other so that thermodynamic properties of SCL contribute to the bending rigidity. Thermodynamic 
properties of SCL are similar to those of Wigner crystal and are anomalous in the sense that the 
pressure, compressibility and screening radius of SCL are negative. This brings about substantial 
negative correction to the bending rigidity. For the case of DNA this effect qualitatively agrees with 
experiment. 

PACS numbers: 77.84.Jd, 61.20.Qg, 61.25Hq 
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I. INTRODUCTION 

Many polymers and membranes are strongly charged 
in a water solution. Among them are biopolymers such as 
lipid membranes, DNA, actin and other proteins as well 
as numerous synthetic polyelectrolytes. In this paper, 
we concentrate on bending of membranes and cylindrical 
polyelectrolytes with fixed uniform distribution of charge 
at their surfaces. For a flat symmetrical membrane, the 
curvature free energy per unit area can be expressed in 
terms of the curvatures c\ and ci along two orthogonal 
axes astl 
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where k is the bending rigidity, kq is the Gaussian rigid- 
ity and S is the membrane surface area. For cylindrical 
and spherical deformations with the radius of curvature 
R c (see Fig. 1) 
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respectively. In general, k = Kq + K e i, where kq is the 
"bare" bending rigidity related to short range forces and 
K e i is electrostatic contribution which is determined by 
the magnitude of surface charge density and the condi- 
tion of its screening by small ions of the water solution. 
Similarly, for a rod-like polymer, such as double helix 
DNA, the change in free energy per unit length due to 
bending is given by 
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where C is the length of the rod, Q = Qo + Qei is the 
bending constant of the rod, which consist of a "bare" 
component, Qo, and an electrostatic contribution Q e j. 
In the worm model of a linear polymer, the persistence 
length, L, of the polymer is related to Q: 
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where Lq is the bare persistent length and L e \ is an elec- 
trostatic contribution to it. In the absence of screen- 
ing, repulsion of like charges of a membrane or a poly- 
electrolyte leads to infinite K e i and L e i. Only screening 
makes them finite. When the surface charge density is 
small enough Debye-Hiickel (DH) approximation can be 
used. For a membrane with the surface charge density 
—a on each side, n e i was calculated au when DH screen- 
ing length r s is larger than membrane thickness h: 
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Here D is dielectric constant of water. 

For cylindrical polyelectrolyte with diameter d much 
smaller than r s , calculations in the DH limit lead to the 
well known Odijk-Skolnick-Fixman formulaQ for the per- 
sistence length: 
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where —rj — irad is the charge per unit length of the poly- 
mer. Eqs. ([3]) and (0) show that, in DH approximation, 
K e i and L e i vanish at r s — so that one can measure Kq 
and Lq in the limit of high concentration of monovalent 
salt. At at r s > 0, the quantities n e i and L e i are always 
positive and grow with r s . However, in many practi- 
cal situations, polyelectrolytes are so strongly charged 
that DH approximation does not work and the nonlinear 
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Poisson-Boltzmann (PB) equation was used to calculate 
K e i and L e i. If counterions Jjiave charge Ze 7 PB equation 
gives, for a thin membraneS 
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where I = Z 2 e 2 / Dk B T is the Bjerrum length with charge 
Z. Eqs. (||), (0), (|), and @ give positive K e i and L e ; in 
agreement with the common expectations that electro- 
static effects can only increase bending rigidity. 

This paper deals with the case of a strongly charged 
membrane or polyelectrolyte with a uniform distribu- 
tion of immobile charge on its surface. It was shown in 
Ref. jsl-fbi] that screening of such surface by multivalent 
counterions with charge Z > 2 can not be described by 
PB equation. Due to strong lateral Coulomb repulsion, 
counterions condensed on the surface form strongly cor- 
related two-dimensional liquid (SCL). Their correlations 
are so strong that a simple picture of the two-dimensional 
Wigner crystal (WC) of counterions on a background of 
uniform surface charge is a good approximation for cal- 
culation of the free energy of the SCL. The concept of 
SCL was used to demonstrate that two charged surfaces 
in the presence of multivalent, counterions attract each 
other at small distancesEIEjO. It was also shown that 
cohesive energy of SCL leads to much stronger counterion 
attraction to the surface than in conventional solutions 
of Poisson-Boltzmann equation, so that surface charge is 
almost totally compensated by the SCLEj. 

In this paper we calculate effect of SCL at the surface 
of a membrane or a polyelectrolyte on its bending rigid- 
ity. When a membrane or polyelectrolyte is bent, the 
density of its SCL follows the changes in the density of 
the surface charge, increasing on one side and decreasing 
on the opposite side of (see fig. 1). As a result the bend- 
ing rigidities can be expressed through thermodynamic 
properties of the SCL, namely two-dimensional pressure 
and compressibility. For two-dimensional one component 
plasma (on uniform background) these quantities were 
found bv-Mante-Carlo simulation and other numerical 
methodslialiZl as functions of temperature. The inverse 
dimensionless temperature of SCL is usually written as 
the ratio of the average Coulomb interaction between ions 
to the thermal kinetic energy kpT 
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where n = a jZe is concentration of SCL. (For e.g., for 
Z = 3 and a = 1.0 e/nm" 2 , T = 6.3). We will show that 
in the range of our interest 3 < T < 15 the free energy, 
the pressure and the compressibility and, therefore, elec- 
trostatic bending rigidities differ only by 20% from those 



in the low temperature limit r — > oo, when SCL freezes 
into WC. General results are given in Sec. III. Here we 
present very simple results obtained in the WC limit: 
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Here a = (2Ze/ v^cr) 1 ^ 2 is the lattice constant of the tri- 
angular close packed WC. The membrane and the cylin- 
der are assumed to be reasonably thick, 2nh ^> a and 
ird ^> a. In contrast with results for DH and PB approxi- 
mations, Kwc an d Lyyc ar e negative, so that multivalent 
counterions make a membrane or a polyelectrolyte more 
flexible. For a membrane with a = 1.0 e/nm -2 , h = 4 
nm at Z = 3 we find that a = 1.7 nm, inequality 2nh 3> a 
is fulfilled and Eq. ( pd| ) yields kwc = — 14fc_eT (at room 
temperature). This value should be compared with typ- 
ical Ko ~ 20 — lOOkpT. For a cylindrical polyelectrolyte 
with parameters of the double helix DNA, d — 2 nm and 
t] — 5.9 e/nm, inequality Ttd ^> a is valid and we ob- 
tain Lwc = —4.9 nm, which is much smaller than the 
bare persistence length Lq = 50 nm. We should, how- 
ever, note that our estimates are based on the use of the 
bulk dielectric constant of water D = 80. For the lateral 
interactions of counterions near the surface of organic 
material with low dielectric constant, the effective D can 
be substantially smaller. (In macroscopic approach it is 
close to D/2). As a result, absolute values of kwc an d 
Ly/c can grow significantly. 

Negative electrostatic contributions to the beiidiag 
rigidity were also predicted in two recent papersEJEM. 
The authors considered this problem in the high temper- 
ature limit where attraction between different points of a 
membrane or a polyelectrolyte is a result of correlations 
of thermal fluctuations of screening atmosphere at these 
points. Such theories describe negative contribution to 
rigidity for Z = 1 or for larger Z but with weakly charged 
surfaces where T < 1. On the other hand, at Z > 3 and 
large ct, one deals with low temperature situation when 
r > 1. In this case the main terms of the electrostatic 
contribution to the bending rigidity are given by Eq. ( |Tl| ) 
and Eq. (p^), which are based on static spatial correla- 
tions of ions. 

We would like to emphasize that, contrary to Ref. [h], 
this paper deals only with small deformations of a mem- 
brane or a polyelectrolyte. We are not talking about a 
global instability of a membrane or polyelectrolyte due 
to self-attraction, where, for example, a membrane rolls 
itself into a cylinder or a polyelectrolyte, as in the case of 
DNA, rolls into a toroidal particleEl Global instabilities 
can happen even when total local bending rigidities are 
still positive. To prevent these instabilities in experiment 
one can wotk with a small area membrane or short poly- 
electrolyteEJ or keep their total bend small by an external 
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force, for example, with optical tweezergPJ'cS. 

It is known that, in a monovalent salt, DNA has a per- 
sistence length L > 50 nm which saturates at 50 nm at 
large concentration of salt. Thus it is natural to assume 
that the bare persistence length La = 50 nm. However, 
it was found in Ref. |2(]-[2^ that a relatively small con- 
centration of counterions with Z — 2,3,4 leads to an 
even smaller persistence length, which can be as low as 
L = 25 — 30 nm. We emphasize that a strong effect was 
observed for multivalent counterions which are known to 
bind to DNA due to the non-specific electrostatic force. 

These experimental data can be interpreted as a result 
of replacement of monovalent counterions with multiva- 
lent ones which create SCL at the DNA surface. As we 
stated before, multivalent counterions should produce a 
negative correction to Lq, although the above calculated 
correction to persistence length is smaller than the ex- 
perimental one. 

This paper is organized as follows. In Sec. II we discuss 
thermodynamic properties of SCL and WC as functions 
of its density and temperature. In sec. Ill and IV we use 
expressions for their pressure and compressibility to cal- 
culate kscl and Lscl and their asymptotic expressions 
Kwc and Lwc- In Sec. V we calculate contributions 
of the tail of screening atmosphere to K e i and L e i and 
show that for Z > 2 and strongly charged membranes 
and polyelectrolytes, tail contributions to the bending 
rigidity are small in comparison with that of SCL. 



II. STRONGLY CORRELATED LIQUID AND 
WIGNER CRYSTAL 

Let us consider a flat surface uniformly charged with 
surface density — a and covered by concentration n — 
o/Ze of counterions with charge Ze. It is well known 
that the minimum of Coulomb energy of counterion re- 
pulsion and their attraction to the background is pro- 
vided by a triangular close packed WC of counterions. 
Let us write energy per unit surface area of WC as 
E = ne(n) where s(n) is the energy per ion. One can 
estimate e(n) as the interaction energy of an ion with 
its Wigner-Seitz cell of background charge (a hexagon of 
the background with charge —Ze). This estimate gives 
e(n) ~ — Z 2 e 2 /Da^ —Z 2 e 2 n 1 / 2 /D. More accurate ex- 
isB 



pression for e(n) 
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where a = 1.96. At room temperature, Eq. ( |13[ ) can be 
rewritten as 
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where a/e is measured in units of nm 2 . 

At a = 1.0 e/nm -2 , Eq. (jlj) gives |e(n)| ~ 7k B T or 
r = 6.3 at Z = 3, and \e(n)\ ~ 13k B T or T = 12 at 
Z = 4. Thus for multivalent ions at room temperature 



we are dealing with the low temperature regime. How- 
ever, it is knownO that due to a very small shear mod- 
ulus, WC melts at even lower temperature: V ~ 130. 
Nevertheless, the disappearance of long range order pro- 
duces only a small effect on thermodynamic properties. 
They are determined by the short range order which 
does not change, .sigiiificantly in the range of our inter- 
est 5 < r < 15M'EJtjO. This can be seen from numeri- 
cal calculationscJ lL3 of thermodynamic properties of clas- 
sical two-dimensional SCL of Coulomb particles on the 
neutralizing background. In the range 0.5 < T < 50, the 
internal energy of SCL per counterion, e(n, T), was fitted 
by 

e(n, T) = fc B T(-l.ir + 0.58r 1/4 + 0.74), (15) 

with an error less than 2% @. The first term on the right 
side of Eq. (|l5|) is identical to Eq. (|l^) and dominates at 
large T. All other thermodynamic functions can be ob- 
tained from Eq. (|l5|). In the next section we show that 
K e i and L e i are proportional to the inverse isothermal 
compressibility of SCL at a given number of ions N 
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where 



P 



-(dF/dS) T = (ne(n,T)+nk B T)/2 
= nk B T(-0.55T + 0.27r 1/4 + 0.87) (17) 

is the two-dimensional pressure, F is the free energy of 
SCL and S = N/n is its area. Using Eq. ( |l7j ) and relation 
dr/dn = r/2n, one finds 

X" 1 = nk B T(-0MT + 0.33r 1/4 + 0.87), (18) 

where the first term on the right side follows from 
Eq. (|l^) and describes WC limit. The last two terms 
give 33% correction to the WC term at T = 5 and only 
12% correction at T — 15. So one can use zero tempera- 
ture, Eq. (|l3|), as first approximation to calculate n e i and 
L e i. This is how we obtained Eq. ( |TT| ) and Eq. ([l2|). 

Eqs. (0) and (||) show that, in contrast with most of 
liquids and solids, SCL and WC have negative pressure 
P and compressibility x- We will see below that anoma- 
lous behavior is the reason for anomalous negative rigid- 
ity K e i and persistence length L e i and positive Gaussian 
rigidity kg, el- The curious negative sign of compressibil- 
ity of two-dimensional electron SCL and WC was first 
predicted in Ref. ^4]. Later it was discovered in magneto- 
capacitance experiments in MOSFETs and semiconduc- 
tor heterojunctionsCJLJ. 

According to Eq. @ = at T = 1.48, P = at 
r = 2.18 and they become positive at smaller T. As one 
can see from Eqs. (14) and (|l0|), at cr ~ 1.0 e/nmT 2 such 
small values of T correspond to Z — 1 . Thus surface layer 
of monovalent ions do not produce large negative K e i and 
L e i in comparison with multivalent ions. For them con- 
ventional results of Eqs. (||), (^), (||), and (^|) related with 
counterions in the long distance tail of screening atmo- 
sphere work better. We will return to this question in 
Sec. V where we discuss the role of these tails. 
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III. MEMBRANE 

We will consider a "thick" membrane for which one can 
neglect the effects of the correlation of SCL on two sur- 
faces of the membrane. If we approximate SCL by WC, 
the energy of such correlations between two surfaces of 
the membrane decay as exp(— 2nh/a), so the condition 
of "thickness" , h ^S> 2ttcl, is actually easily satisfied for a 
strongly charged membrane. 

Let us first write the free energy of each surface of the 
membrane as 



F = Nf(n,T) 
where f(n,T) is the free energy per ion. 



(19) 
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a) 



b) 



FIG. 1. Bending of membrane (the curvature has been ex- 
aggerated). For simplicity, the WC case is depicted, a) A 
thick membrane. The right WC is compressed while the left 
WC is stretched. For thick membranes, this is the dominant 
cause of the change in free energy, b) A very thin membrane. 
Only one Wigner-Seitz cell is shown. Due to finite curvature 
of the surface, the distance from any point of the Wigner-Seitz 
cell to the central ion is shorter than that in the flat config- 
uration. For thin membranes, this is the dominant cause of 
free energy change. 

When a membrane is bent (see Fig. la), the surface 
charge on the right side is compressed to a new den- 
sity n R > n, while the surface charge on the left side is 
stretched to nj, < n. Since the total charge on each sur- 
face is conserved, this change in density leads to a change 
in the free energy of each surface: 
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in which we kept only terms up to second order in 
Sn L:R = n L:R - n. 

Using the definitions ( |l7| ) and ([l6|) for the pressure and 
the compressibility of 2D systems 
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Eq. (PQ ) can be rewritten as 
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So, the total change in the free energy of the membrane 
per unit area is 
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In the case of cylindrical geometry, keeping only terms 
up to second order in the curvature R ~ , we have 
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Substituting Eq. @ into Eq. (||), we get 

-h 2 R~ 2 . 



(25) 



S 4 X 

Similarly, in the case of spherical geometry we have 
R c 



(26) 



n L ,R 



R c ± h/2 



h 3h 2 \ 

l *K + im) n (27) 



and 



S~ 



1 P 

X ~ 2" 



h 2 R- 



(28) 



Comparing Eq. @) and (||) with Eq. (§) and (|), we 
obtain general expressions for the electrostatic contribu- 
tion to the bending rigidity 
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For example, in the case of low surface charge density, 
DH approximation can be used to geto 



/(n,T) = 27r— n-V, , (30) 

from which, we can easily get a generalization of Eq. (|^) 
for a "thick" membrane (h 3> r s ) 



k dh = 2n—h 2 r s , Kg,DH 
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(31) 



In the case of high surface charge density we study in 
this paper, a SCL of multivalent counterions resides on 
each surface of the membrane. The expressions for the 
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pressure and the compressibility given by Eqs. ([L7|) and 
( |l8| ) can be used to calculate the bending rigidity: 



kscl = —k B T(-0.83T + 0.33T^ 4 + 0.87) , (32) 
k g ,scl = -^fc s T 1 (-0.55r + 0.27r 1 / 4 + 0.87) . (33) 



In the limit of a strongly charged surface (r> 1), the 
first term in Eqs. (|32| ) and ([53]) dominates, the free en- 
ergy of SCL is close to that of WC. Using Eq. ( |io| ) one 
arrives at Eq. (11) for the bending rigidity in the WC 
limit. 

As already stated in Sec. 1, for T > 3, Eqs. (32), 
( |33| ) give a negative value for the bending modulus and 
a positive value for the Gaussian bending modulus. In 



other words, multivalent counterions make the membrane 
more flexible. This conclusion is opposite to the standard 
results obtained by mean field theories (Eqs. (g), (§), 
(|3l|)) where electrostatic effects are known to enhance the 
bending rigidity of membranes (n e i > and nc,ei < 0). 
Obviously, this anomaly is related to the strong correla- 
tion between multivalent counterions condensed on the 
surface of the membrane, which was neglected in mean 
field theories. 

We can also look at Eqs. (|3|) and ( |ll| ) from another 
interesting perspective: apart from a numerical factor, 
Eq. (31) is identical to Eq. ( |TT| ) if we replace r s by —a. 
So the WC of counterions has effect on bending proper- 
ties of the membrane as if one replaces the normal 3D 
screening length of counterions gas by a negative screen- 
ing length of the order of lattice constant. Such negative 
screening length of WC or SCL has been derived for the 
first time in Ref. |27j. It follows from the negative com- 
pressibility predicted in Ref. Q and observed in Refs. |25| 
and |§ 

Until now we have ignored the effects related to Pois- 
son's ratio ap of the membrane material. We are talking 
about the bending induced increase of the thickness of 
the compressed (right) half of the membrane, simultane- 
ous decrease of the thickness of its stretched (left) half, 
and the corresponding shift of the neutral plane of the 
membrane (the plane which by definition does not expe- 
rience any compression or stretching) to the left from the 
central plane. These deformations can be found following 
Ref. ^8] and lead to additional term aph 2 / (1 — ap)R 2 in 
the right side of Eq. (p5h. It gives for the bending rigidity 
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So, for example, at ap = 1/3, the second term of Eq. ( |34| ) 
gives a 33% correction to Eq. (|TT|) . 

According to Eqs. ©, (g^T^l) « eJ = at h = 0. 
This happens because in this limit two SCL merge into 
one, whose surface charge density remains unchanged af- 
ter bending. Nevertheless, there is another effect directly 
related to the curvature of SCL. It can be explained by 



concentrating on one curved Wigner-Seitz cell (see Fig. 
lb). One can see, that due to the curvature, points of 
the background come closer to the central counterion of 
the cell in the three-dimensional space where Coulomb 
interaction operates. As a result, the energy of SCL goes 
down. In the Wigner-Seitz approximation, where energy 
per ion of WC is approximated by its interaction with the 
Wigner-Seitz cell of the background charge, we obtain 
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We see that this effect also gives anomalous signs for elec- 
trostatic contribution to rigidity in the WC limit, but 
with a very small numerical coefficient. Also note that, 
as in the thick membrane case, we can obtain Eq. ( ^H|) for 
a thin membrane by replacing r s in Eq. (g) by a negative 
screening radius of WC with absolute value of the order 
n. 



IV. CYLINDRICAL POLYELECTROLYTES 



In this section, we study bending properties of cylin- 
drical polyelectrolytes with diameter d and linear charge 
density rj (see Fig. 2). As in the membrane problem, 
we will assume that the cylinder is thick, i.e. its cir- 
cumference nd is much larger than the average distance 
a between counterions on it surface. The calculation is 
carried out exactly in the same way as in the case of thick 
membrane. The only difference is that, instead of sum- 
ming the free energy of two surfaces of the membrane, 
we average over the circumference of the cylinder. 

Let us denote by the local density at an angle <f> on 
the circumference on the cylinder (see Fig. 2a). Before 
bending = n = n/ndZe, after bending it changes to a 
new value 
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Using Eq. (gj) the free energy per unit length of the 
polymer can be written as 
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where we keep terms up to second order in the curvature 
R„ . 
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Obviously, due to the expansion in y direction, the 
correction to energy is not as strong as in the membrane 
case. For example, at <jp — 1/3, Eq. ( f4l| ) gives only 3% 
correction to Eq. ( [14). 

According to Eqs~(||) and (|l|), at d — 0, K e i vanishes. 
In this limit, we have to directly include the curvature 
effect on one dimensional SCL as shown in Fig. 2b. As 
already mentioned in the previous section, after bending, 
points on a Wigner-Seitz cell come closer to the central 
ion, which lowers the energy of the system. This effect 
can be calculated easily in the WC limit. Let's consider 
the electron at the origin, its energy can be written as 



FIG. 2. Bending of cylindrical polyelectrolytes. a) A thick 
cylinder. Rigidity is mostly determined by the change in den- 
sity of SCL. b) A thin cylinder. The curvature effect, is the 
dominant cause of change in free energy. 



Comparing Eq. ( j37| ) with Eq. (||), (||), one can easily 
calculate the electrostatic contribution to the persistence 
length 
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In the case of highly charged polymer, a SCL of counteri- 
ons resides on the polymer surface. For a thick cylinder, 
the SCL is locally flat and we can use the numerical ex- 
pression (|l8| ) for \~ x to obtain 



Lscl = ^nd 3 (-0.83r + 0.33r 1/4 + 0.87) 



(39) 



Again, we see that correlations between counterions on 
the surface of a polymer lead to a negative electric con- 
tribution to persistence length for T > 1.5. In the WC 
limit r ^ L the first term in Eq. ( |39"| ) dominates, and 
using Eq. ( [Toj ) one can easily obtain Eq. (|l2]). 

As in the membrane case, for simplicity, in writing 
down Eqs. ([56|), we have ignored the effect of finite value 
of the Poisson's ratio of the polymer material. In mem- 
branes, this effect result in a gain in energy due to the 
shift of the neutral plane toward the convex (stretched) 
sides. For a cylinder, there is an additional expansion in 
the y direction (Fig. 2) which reduces the change in sur- 
face charge density, hence compensates the above gain. 
These deformations can be found following Ref. |2^ and 
lead to a correction to Eqs. (pq) 
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This gives, for the persistence length, 
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where ri = ia and s is the contour distance from our ion 
to an lattice point i and the element ds of the background 
charge. In the straight rod configuration the space dis- 
tant is the same as the contour distance, however after 
bending they change to 

r^r,(l-r 2 /24i? c 2 ) , s' ~ s(l - s 2 /24R 2 c ) . (43) 

Using these new distances to calculate the energy of the 
bent rod and subtract Eq. (|4|) from it, one can easily 
calculate the change in energy due to curvature and the 
corresponding contribution to persistence length: 



j thin 



I 

96 



which is negative and very small. For e.g., for Z 
-0.065 nm and —0.116 nm respectively. 



(44) 



3, 4, 



t thin 



V. CONTRIBUTIONS OF THE TAIL OF THE 
SCREENING ATMOSPHERE 

In previous sections, we calculated the contribution of 
a SCL of counterions condensed on the surface of a mem- 
brane or polyelectrolyte to their bending rigidity. We as- 
sumed that charge density a is totally compensated by 
the concentration n = a / Ze. Actually, for example, for a 
membrane, some concentration, N(x), of counterions is 
distributed at a distance x from the surface in the bulk of 
solution (we call it the tail of the screening atmosphere) . 

The standard solution of PB equation for concentra- 
tion N(x) at N(oo) = has a form 



N{x) 



1 



1 



2ttZ (A + x) 2 



(45) 



where A = Ze/(2irla) is Gouy-Chapman length. At 
r 3> 1, correlations in SCL provide additional strong 
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binding for cotrnterions, which dramatically change the 
form of N(x)t3. It decays exponentially at A <C x <C Z/4, 
and at x ^> Z/4 it behaves as 



AT(a;) 



2ttZ (A + rr) 2 



(46) 



Here A = Ze/(2Trla*) is an exponentially large length 
and a* is the exponentially small uncompensated surface 
charge density at the distance ~ Z/4. In any realistic sit- 
uation when N(oo) is finite or a monovalent salt is added 
to the solution, Eqs. ( fl5| ) and (^) should be truncated 
at the screening radius pp . Then the solution of the stan- 
dard PB equation givesa Eq. (^) at r s > A or Eq. @ at 
r s -C A. In the case of SCL, for realistic values of r s in 
the range 1/4 < r 5 « A, we obtain a contribution of the 
tail similar to Eq. (||) 



n t = 3tt- 



D 



(47) 



At reasonable values of r s , this expression is much smaller 
than Kwc due to very small values of the ratio a* /a. 

Now we switch to a cylindrical polyelectrolytfL In this 
case, the solution of the PB equation is knowncE to con- 
firm the main features of the Onsager-ManningE2l picture 
of the counterion condensation. This solution depends 
on relation between \rj\ and rj c = Ze/l. In the case inter- 
esting for us, 1 77 1 3> 77c, the counterion charge \rf\ — rj c is 
localized at the cylinder surface, while the charge ?7 C , is 
spread in the bulk of the solution. This means that at 
large distances the apparent charge density of the cylin- 
der, r] ai equals — i] c and does not depend on 77. Eq. (j^) 
can actually be obtained from Eq. (0) by substituting r\ c 
for 77. 

It is shown in Ref . |lj that at T ^> 1 , the existence of 
SCL at the surface of the cylinder leads to substantial 
corrections to the Onsager-Manning theory. Due to ad- 
ditional binding of counterions by SCL |?7 a | < |?7 C | and is 
given by the expression 



ln[iV(0)/iV(oo)] 
: ln(4r s /Z) 



(48) 



where N(0) is exponentially small concentration at the 
distance r > Z/4 from the cylinder axis, used in Ref. |l4|as 
a boundary condition for PB equation at x = 0. There- 
fore, one can obtain for the tail contribution, the estimate 
from the above using Eq. (||). For Z = 3 and r s = 5 nm 
this gives L t < 1 nm. For DNA, this contribution is much 
smaller than L$cl — — 5 nm. 



VI. CONCLUSION 

We would like to conclude with the discussion of ap- 
proximations used in this study. First, we assumed that 
the surface charges are immobile. This is true for rigid 
polyclcctrolytcs, such as double helical DNA or actin, as 



well as for frozen or tethered membranes. But if the 
membrane is fluid, its charged polar heads can move 
along the surface. In this case surface charges can ac- 
cumulate near a Z-valent counterion and screen it. Such 
screening creates short dipoles oriented perpendicular to 
the surface. Interaction energy between these dipoles is 
much weaker than the correlation energy of SCL. There- 
fore it produces negligible contribution to the membrane 
rigidity. The mobility of the charged polar heads elimi- 
nates effects of counterion correlation only in the situa- 
tion where the membrane has polar heads of two differ- 
ent charges, for example, neutral and negative ones. In 
such a membrane, the local surface charge density can 
grow due to the increase of local concentration of nega- 
tive heads. But if all of the closely packed polar heads 
are equally charged their motion does not lead to redis- 
tribution of the surface charges. Then our theory is valid 
again. 

Another approximation which we used is that the sur- 
face charge is uniformly smeared. This can not be exactly 
true because localized charges are always discrete. Nev- 
ertheless our approximation makes sense if the surface 
charges are distributed evenly, and their absolute value 
is much smaller than the counterion charge. For example, 
when the surface charged heads have charge -e and the 
counterion charge is Ze 3> e, then the repulsion between 
counterions is much stronger than their pinning by the 
surface charges. At Z > 3 we seem to be close to this 
picture. On the other hand, if the surface charges were 
clustered, for example, they form compact triplets, the 
trivalent counterion would simply neutralize such clus- 
ter, creating a small dipole. Obviously our theory would 
over-estimate electrostatic contribution to the bending 
rigidity in this case. 

All calculations in this paper were done for point like 
counterions. Actually counterions have a finite size and 
one can wonder how this affects our results. Our results, 
of course, make sense only if the counterion diameter is 
smaller than the average distance between them in SCL. 
For a typical surface charge density, a = 1.0 e/nm~ 2 , the 
average distance between trivalent ions is 1.7 nm, so that 
this condition is easily satisfied. The most important cor- 
rection to the energy is related to the fact that due to 
ion's finite size, the plane of the center of the counterion 
charge can be located at some distance from the plane of 
location of the surface charge. This creates an additional 
planar capacitor at each surface and results in a positive 
contribution to the bending rigidity similar to Eq. (|3l| ) 
which can compensate our negative contribution. On the 
other hand, if the negative ions stick out of the surface 
and the centers of counterions are in the same plane with 
centers of negative charge this effect disappears. 

In general case, one can look at this problem from an- 
other angle. Let us assume that the bare quantities kq 
and Lq are constructively defined as experimental values 
obtained in the limit of a high concentration of monova- 
lent counterions. Let us also assume that the distances of 
closest approach of monovalent and Z-valent counterions 
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to the surface are the same. This means that the pla- 
nar capacitor effect discussed above is already included 
in the bare quantities kq and Lq. Then the replacement 
of monovalent counterions by Z-valent will always lead 
to Eq. @ and Eq. @. 

In summary, we have shown that condensation of mul- 
tivalent counterions on the surface of a charged mem- 
brane or polyelectrolytc happens in the form of a strongly 
correlated Coulomb liquid, which closely resembles a 
Wigner crystal. Anomalous properties of this liquid lead 
to the observable decrease of the bending rigidity of a 
membrane and polyelectrolyte. 
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